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Spin-dependent transport through an interacting single-level quantum dot coupled to ferromag- 
netic leads with non-coUinear magnetizations is analyzed theoretically. The transport properties 
and average spin of the dot are investigated within the nonequilibrium Green function technique 
based on the equation of motion in the Hartree-Fock approximation. Numerical results show that 
Coulomb correlations on the dot and strong spin polarization of the leads significantly enhance 
precession of the average dot spin around the effective molecular field created by the external elec- 
trodes. Moreover, they also show that spin precession may lead to negative differential conductance 
in the voltage range between the two relevant threshold voltages. Nonmonotonous angular variation 
of electric current and change in sign of the tunnel magnetoresistance are also found. It is also 
shown that the diode-like behavior in asymmetrical junctions with one electrode being half-metallic 
is significantly reduced in noncoUinear configurations. 

PACS numbers: 72.25.Mk; 73.63.Kv; 73.23.Hk 
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I. INTRODUCTION 

Current interest in electronic transport and spin ef- 
fects in mesoscopic tunnel junctions is stimulated by their 
possible applications in microelectronics and spintronics 
devices |^. One of the most widely studied spin- 
dependent effects in magnetic tunnel junctions is the tun- 
nel magnetoresistance (TMR) . This phenomenon appears 
as a change in the junction resistance when magnetic 
moments of external electrodes rotate from parallel align- 
ment to non-collinear one (or to antiparallel alignment in 
a particular case) . Such a rotation of magnetic moments 
may be induced, for instance, by an external magnetic 
field. The TMR effect may occur in planar junctions 
HiQiUli mesoscopic double-barrier junctions [sSHjIS' 
granular systems [Tol and others. When the cen- 
tral part of a double-barrier junction is sufficiently small, 
the interplay of discrete charging by single electrons and 
spin dependence of tunnelling processes can lead to addi- 
tional interesting features in the corresponding transport 
characteristics jS,, A »2lJ • 

Up to now, most of theoretical works on transport 
through quantum dots (QDs) coupled to ferromagnetic 
electrodes was limited to collinear, i.e., parallel and an- 
tiparallel magnetic configurations 0, 0, 0, 0, 0, 
Il8| . It is only very recently when transport in sys- 
tems with non-collinear magnetizations was addressed 
[l^ I20II21II23 I . In particular, it has been shown that the 
diode-like features in transport characteristics of systems 
with one electrode being half-metallic are significantly re- 
duced when magnetic moments of the electrodes become 
non-collinear • However, this behavior was studied 
only in the sequential tunneling regime. 

In recent papers spin precession in electron 

tunnelling through an interacting quantum dot was stud- 



ied theoretically in the first order approximation with 
respect to the tunneling Hamiltonian. Such a preces- 
sion takes place when magnetic moments of the leads are 
non-collinear, and occurs because an electron entering 
the dot in a tunnelling event is subject to an effective ex- 
change field, which exerts a torque on the electron spin 
and makes the spin precesses by a certain angle before 
the electron leaves the dot. In the first order approx- 
imation, the spin precession is driven by the Coulomb 
interaction on the dot (described by the Hubbard cor- 
relation parameter U) and disappears when U tends to 
zero. 

In this paper we extend the earlier descriptions pol 
I21I of electron tunneling through QD with non- 
collinear magnetizations by going beyond the first order 
approximation. In order to calculate the tunneling cur- 
rent and spin precession, we employ the non-equilibrium 
Green function technique and limit considerations to the 
Hartree-Fock approximation. Hence, besides the sequen- 
tial tunneling also the contribution due to higher order 
tunneling processes is included in the description. Nu- 
merical results show, that spin precession also exist in the 
limit of U — 0. However, Coulomb correlations signifi- 
cantly enhance the precession. This, in turn, may lead to 
negative differential conductance in a certain bias voltage 
range. We predict that symmetry of the junction, bar- 
rier height, and spin polarization of magnetic leads may 
significantly influence spin dependent characteristics of 
the dot. In particular, it is shown that the diode- like be- 
havior in asymmetrical junctions with one half-metallic 
electrode is partially suppressed in non-collinear config- 
urations. However, the suppression is much less evident 
than that obtained within a simplified theory neglec ting 
the exchange interactions between the dot and leads ■ 

The paper is organized as follows. In Section II we 



2 



describe model of the system. Theoretical method is 
described in Section III, where the equation of motion 
method is used to derive nonequilibrium Green functions 
of the dot. Transport characteristics are calculated in 
Section IV, whereas relevant numerical results on tun- 
nelling current, magnetoresistance, and spin precession 
are presented and discussed in Section V. Finally, sum- 
mary and general conclusions are in Section VI. 



II. MODEL 

We consider a single-level QD coupled to two ferro- 
magnetic metallic leads by tunnelling barriers. Magnetic 
moments in the external leads lie in a common plane and 
form an arbitrary angle ^p. To describe electron spin wc 
will use the local and global quantization axes. The lo- 
cal axes are determined by the local spin polarization in 
the leads. The global quantization axis (the axis z in 
our case) is assumed to coincide with the local one in 
the left electrode. Spin projection on the local quantiza- 
tion axes will be denoted as /3 = + for majority electrons 
and (3 = — for minority electrons, whereas projection 
on the global quantization axis will be denoted as a =t 
and a —[. Axis y of the coordinate system is normal to 
the plane xz determined by the spin polarizations of the 
leads. Geometry of the device and the orientation of the 
coordinate system are shown schematically in Fig.l. 
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FIG. 1: Schematics of the system considered in this paper. 
The coordinate systems used to describe states of the dot is 
also shown. 

The whole system can be described by Hamiltonian of 
the general form 



H — Hi + Hr + Hd 



(1) 



The term H^, describes the left (i^ = I) and right {v — r) 
electrodes in the non-interacting quasi-particle approxi- 
mation, 

■f^-^ = X! X! ^wo-thpo-^kf}, (2) 

k 13=+ - 

where is the single-electron energy in the v-fh. elec- 
trode for the wavevector k and spin /3, whereas a^j,^ and 
o-vkp are the corresponding creation and annihilation op- 
erators. The term Hd in Eq.(l) describes QD, 



where = c+Co- is the occupation operator, Ed denotes 
the energy of the discrete level, U is the electron corre- 
lation parameter, whereas and Co- are the correspond- 
ing creation and annihilation operators for electrons with 
spin orientation a {[). Finally, the tunnelling term, 
Ht, in Eq.(l) takes the form 



Ht = Hi + m 



(4) 



where the first term describes tunnelling through the left 
barrier. 

Hi = ^(T^+a+ + CT + T^a+.cx) + h.c, (5) 

k 

whereas the second term corresponds to tunnelling 
through the right barrier, 

HI - y{m+a+k+ cos(^/2) - n_a+^_ sin(^/2)]cT 



fc 



+ K-<fe- cos(^/2) + T,r+a+ ^ sin((p/2)]cx} + h.c, (6) 

and h.c. stands for the hcrmitian conjugate terms. The 
tunneling terms Hi and HI have different forms because 
the corresponding local and global axes are parallel for 
the left electrode and noncoUinear for the right one. 



III. GREEN FUNCTIONS OF THE DOT 

To calculate electric current in nonequilibrium situ- 
ations we will make use of the nonequilibrium Green 
function defined on the Keldysh contour [2^. The 
causual Green function of the dot is defined as Gaa' (e) = 
{{ca\c^i))<L- Writing equation of motion for {{ccr\c^i))<i, 
one arrives at 

(e-ed)((c^|c+)), = 5„a' 

+ Y}Ttp{M c+ )). + n}>{{arkp\c+,)), cos(^/2) 

k 

-/?T*!:^((a.fc_^|c+)),sin((^/2)] + C/((c.7w|c+))„ 

(7) 

where (3 — + for u =] and [3 = — for <j =[. Applying 
equation of motion to the four new Green functions on 
the r.h.s. of Eq.(7), one finds 



(6-4^)((a,fe0|c+)), = r^^((c.|c+)), 



(8) 



(e - el^){{ark0\ c+ )). = T^^((c.| c+ )), cos(<^/2) 
+/3r^^((c.,|c+)),sin(^/2), (9) 

6 - el_p){{ark-0\ c+)). = n_p{{c^„\ c+)),cos(^/2) 
-/3r,'-_^((c.|c+))eSin((^/2), (10) 



Hd = '^edC^Cc +Un^ni, 



(3) 
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{e-ea-U){{c,n.,\c+.)), = {{c,n.,, c+ }) + 5^[T,*^((a;fc;3| c+ )), + T,*^^((a,fc;3 n_<.| c+ )), cos(<^/2) 

k 

-PT^HfAia-rk-pn^al c+,}),sin((^/2) - Tl_i^{{caa+_f^c^a\ c+,))e - T^ipHcaaik-ffct^l c+,))e 
-n_p{{c^a+^_^c^,\ c+ )), cos(^/2) - f3T^p{{c,a+^^c^,\ c+ )), sin(^/2) 

-T*Lfs{{c.ark-0ctJ c+ )), cos(^/2) - /3T*^((c,a,fe^c;j c+ ))e sin(^/2)]. (11) 



Now, the Hartree-Fock decoupling scheme is apphed to where 
the higher-order Green functions generated on the r.h.s. 
of eq. (11), Eio)(g) 



{{a^k±fjn-a\c'^,))e {n-a){{auk±fj\Cy))e, (12) 



{{Caauk±0ct„\c+,)) {ct„Ca){{a^k±0\4>)), (13) 



^k±\ 



k ^ ^k± 

1 2 iT'r 1 2 



(19) 



fc±/3 ' 



(14) 



which closes the set of equations (7)-(ll) and allows to 
find solution for the causual Green functions Gaa-'i^)- 
Here, (...) means the quantum statistical average value 
of the appropriate operator. 

The general solution for the Green function Gcra' (e) in 
the Hartree-Fock approximation may be written in the 
matrix form as 

G(e) = [l-g(e)SW(e)]-ig(e) = [l-gW(e)S(e)]-ig(")(6), 

(15) 

where gl^J,{e) ~ S^^i{e — e^)^^ is the Green function of 
the dot in the absence of both Coulomb interaction and 
coupling to the leads, whereas g(e) is the corresponding 
Green function of the dot decoupled from the electrodes 
but with the Coulomb interaction included. 



5(7(7 (e) = 



e-ed-U{l- {n-a)) 
{e - £d){e - £d - U) 



(16) 



k+\ 



sm<^ . 



(20) 



As follows from Eq.(15), the full self-energy S(e) is de- 
termined by T!i^^\e) and g~^(e) via the relation 



S(e) = S(0)(e)+gW-He)-g-'(e). 



(21) 



The retarded and advanced Green functions G^^^)(e) 

can be found as G^^^^{() = G(e ± 177). Similarly, one 
can find Y,"(^^{e) and H^^^^^^^e). From Eqs (18) to 
(20) follows that the retarded self-energy S(°)^(e) has the 
general form (18) with the components given explicitly 
by the formulae 



(e) 



1^ /D + elW 
TT \D -eVi+e 



+ i 



-- [ri(e)cos^(<^/2) +r^^(e)sin2(v,/2)] 



1 In ( D + ^'^r-i 



(22) 



{e - £d){e - Ed - Uy 



(17) 



with {n_aa) = {ct^Ca). In Eq.(15) S(°)(e) is the self- 
energy of the noninteracting [U = 0) quantum dot, 
whereas S(e) is the full self-energy including contri- 
butions from the coupling to the leads and from the 
Coulomb interaction. The self-energy S(°)(e) takes the 
form. 



S(°)(e) I]^'^(e) 



(0)/ 



(18) 



.(0)7? 



where 



(e) = --[r;(e)-r!l(e)]sin<p 
1 (D + eVr-e 



ri(e)=27r^|T,^±|2,5(e-e^±) 



(23) 



(24) 



for y = l,r. It has been assumed that the lower and 
upper edges of the electron band at zero bias are at —D 
and D, respectively. 
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In the following we assume 

r^(e) = r'± = ro(i±p/) (25) 

and 

Tl{e)^Tl^aTo{l±Pr) (26) 

when e is within the electron band and zero otherwise. 
The parameters p ; and p r describe the spin asymmetry of 
the coupling to the left and right electrodes, respectively, 
Fq is a constant, and the parameter a takes into account 
asymmetry between coupling of the dot to the left and 
right electrodes. 

The key problem is to find the correlation Green func- 
tion G<(e). This can be calculated from the Keldysh 
equation, 

G<(e)=G^(e)S<(e)G^(e), (27) 

provided the lesser self-energy S"^(e) is known. How- 
ever, finding S^(e) for an interacting {U > 0) system is 
a difficult tusk and one has to use some approximate so- 
lutions. One approach is the so-called Ng ansatz, accord- 
ing to which the lesser (greater) self energy S<'^>)(e) is 
proportional to the self-energy of the corresponding non- 
interacting system, S<(>)(e) — AS(°)<'>)(e), where A 
can be be determined on taking into account the relation 
S<(e) - S>(e) = S^(e) - S'^(e). Taking into account 
the above and Eq.(21) one finds A = I, where I is the 
unit matrix. Thus, in the Hartree-Fock approximation 
one finds 

S<(e) = -^[S(°)^(e) - SW-^(e)]/.(e), (28) 

where /^.(e) is the Fermi-Dirac distribution function for 
the i^-th electrode, = 1/{1 -|- exp[(e — ii^)/kBT]}, 

with the electrochemical potentials /z/ = eVi = eV/2 and 
^ir = eVr = —eV/2. (The energy is measured from the 
Fermi level of the leads in equilibrium.) In the Hartree- 
Fock approximation the relation (28) can also be ob- 
tained from the equation of motion method. 

The average values of the occupation numbers (n^) — 
{c^c„) and {ria^^) = (c+c_o.), which enter the expres- 
sions for Green functions, have to be calculated self- 
consistently by using the formulae 

/■+°° de 

K)=Imy^ ^G<Ae) (29) 

and 

/ + 00 7 
^ ^G<,,(e). (30) 



IV. TRANSPORT CHARACTERISTICS 

Having found the Green functions, one can calculate 
electric current, spin accumulation, and spin precession 
on the dot. To do this we calculate the average values of 
all the three spin components, which are related to the 
diagonal and off-diagonal occupation numbers (Eqs (35) 
and (36)). The corresponding relations may be written 
as (spin components are measured in the units of h), 

(5,) = (nT~nx)/2, (31) 



{Sy)=lm{nn), (32) 



(S,) - Re{nn), (33) 

where n-f, n^, and are calculated self-consistently, as 
described above. 

In turn, electric current flowing from the i^-th lead to 
the dot is given by the formula |24| 

"^"^ = f {r.(G<(e) + fA^)[G%) - G%)])} , 

(34) 

with 
and 

Tr=Toa(^+P^.''°''^ ^ P^''""^ ). (36) 

\ Pr Sm ip 1 — Pr cos if J ^ ' 

Thus, taking into account Eqs (22)- (34), together with 
Eqs (40)-(42), one obtains the final symmetrized expres- 
sion for electric current, J = (1/2) (J/ — J,.), in the form 



J=^J^ de[Me)~Me)]j{e), (37) 

where 
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j(e) = 2(l+pO(l+PrCOS(^)Gf^(e)Gf<*(e) + 2(l-pi)(l-PrCOS</.)G«(e)Gf<*(e) 

+ [(1 -Pr-cos^) + (1 +PrCOS^)][Gf^(e)Gf;(e) + Gf^ (e)Gf^* (e)] 

+ (1 +pOP.sin(^{Gf^(e)[Gf;(e) + Gff (e)] + Gff (e)[Gf^(6) + Gf^(e)]} 

+ (1 -pOP.sin^{Gf.(6)[Gfne) + Gff (6)] +Gfne)[Gf,(6) +Gf,(6)]}. (38) 



The current formula derived above can be applied to 
any magnetic configuration of the system, and thus can 
be used to determine TMR. Generally, TMR is described 
quantitatively by the ratio 

TMi.^^M,^^^^^^fc5-^, (39) 



Ri^ = 0) 



J(^) 



where J{ip) is the electric current flowing through the 
system when the angle between spin polarizations of the 
leads is cp, and R{(p) is the corresponding electrical resis- 
tance. 



V. NUMERICAL RESULTS 
A. Symmetrical junctions 

In a symmetrical case both barriers are identical, a — 
1 , and the electrodes are made of the same ferromagnetic 
material, pi = Pr- Consider first electronic transport in a 
symmetrical junction with fully polarized (half-metallic) 
external electrodes, pi — Pr = 1, and with the dot level 
above the Fermi level of the electrodes at equilibrium, 
Sd > 0. 

Bias dependence of electric current and the corre- 
sponding differential conductance are shown in Fig. 2(a) 
and Fig. 2(b) for selected values of the angle (f. The 
current-voltage curve for parallel configuration [ip — 0) 
reveals typical step-like characteristics. Below the first 
(lower) threshold voltage the dot is empty and thus se- 
quential contribution to electric current is exponentially 
suppressed. The first step in the current occurs at the 
bias, where the discrete level Ed crosses the Fermi level of 
the source electrode (the dot can be occupied by a single 
electron), whereas the step at a higher voltage (higher 
threshold) corresponds to the case when Sd + U crosses 
this Fermi level (the dot may be doubly occupied). 

In non-coUinear configurations a monotonous suppres- 
sion of the tunneling current with increasing if takes place 
in the whole bias voltage range, and the current disap- 
pears for if = TT. This is a typical (perfect) spin-valve 
effect. Suppression of electric current is due to an elec- 
tron residing on the dot, whose spin orientation prevents 
it from tunneling to the drain lead. In the extreme case 
of antiparallel configuration, if ~ tt, the electron that 
has tunneled to the dot from the fully polarized source 
electrode blocks transport through the junction since it 



cannot tunnel further to the oppositely polarized drain 
lead. This scenario holds as long as spin-flip relaxation 
processes are absent. 

The steps in current-voltage characteristics give rise to 
the narrow peaks in differential conductance displayed in 
Fig. 2(b), which occur at the lower and higher threshold 
voltages. Apart from this, dependence of electric cur- 
rent on magnetic configuration of the system leads to 
the TMR effect, defined quantitatively by Eq.(39), and 
is shown in Fig. 2(c). The effect increases with increas- 
ing angle ip and tends to infinity when (p —>■ n (therefore 
there is no curve in Fig. 2(c) for ip — tt). 

An interesting feature of the current-voltage charac- 
teristics is the negative differential conductance, which 
may occur in non-coUinear configurations between the 
lower and higher threshold voltages (between the corre- 
sponding two peaks in the differential conductance). The 
negative differential conductance corresponds to some en- 
hancement in TMR, as can be seen in Fig. 2(c). The 
enhancement is particularly significant for rather large 
values of the angle ip (smaller that tt). Physical ori- 
gin of this feature becomes clear, when taking into ac- 
count bias dependence of the average spin on the dot in 
nonequilibrium situation (spin accumulated on the dot). 
As shown in Fig. 3(a), the absolute value of the average 
(Sz) increases with increasing ip and almost vanishes in 
the parallel configuration. In turn, the average value of 
the x-componcnt, (Sx), vanishes for both parallel and an- 
tiparallel alignments and is nonzero for canted configura- 
tions, as shown in Fig. 3(b). Similarly, the average value 
of (Sy) also vanishes in the coUinear configurations and is 
nonzero in the noncoUinear ones. One should point here, 
that neither initial spin state in the source electrode, nor 
the final state in the drain electrode have nonvanishing 
y-component (perpendicular to the plane determined by 
spin polarizations of the two leads). What is then the 
reason of nonvanishing {Sy)7 This can be accounted for 
by taking into account the fact, that an electron resid- 
ing on the dot experiences a certain exchange field re- 
sulting from coupling between the dot and leads, which 
effectively acts as a local magnetic field 22]. Strength 
of this molecular field and its orientation with respect 
to the global quantization axis depend on the applied 
bias voltage and on the angle between the magnetic mo- 
ments of the leads. Thus, if an electron that has tunneled 
from the source lead resides sufficiently long time on the 
dot level, its spin experiences a torque due to the ex- 
change interaction, which results in precession of the av- 
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FIG. 2: Bias dependence of electric current (a), differential 
conductance (b), and tunnel magnetoresistance (c), calcu- 
lated for indicated values of the angle ip. The inset in (a) 
shows electric current between the lower and higher threshold 
voltages for = 27r/3. The parameters assumed for numeri- 
cal calculations are: = 0.1 eV, U = 0.4 eV, To = 0.01 eV, 
Pi = Pr = 1, a = 1 and r = 100 K. 



erage spin around the molecular field, and consequently 
in a nonzero average value of the transverse component 
{Sy). As shown in Fig. 3(c), this precession- induced com- 
ponent is significant in non-coUinear cases (see the case 
of — 27r/3) and in the voltage range between the two 
threshold voltages. Just this enhanced precession is as- 
sociated with a decrease in electric current, which leads 
to negative differential conductance. 

As follows from the above discussion, the magnitude 
of {Sy) is a measure of the spin precession induced by 
the effective exchange field. In Fig. 4(a) we show {Sy) for 
different values of the Coulomb correlation parameter U . 
The curve for U — Q indicates that spin precession in 
noncollinear configurations takes place also when there 
is no Coulomb interaction between electrons on the dot. 
However, as follows from Fig. 4(a), the presence of such an 
interaction enhances the spin precession and also extends 
the voltage range where the precession is significant. 



FIG. 3: Bias dependence of the average spin components; {Sz) 
(a), {Sx) (b) and {Sy) (c) for indicated values of the angle ip. 
The other parameters are as in Fig. 2. 



The spin precession also depends on the spin polariza- 
tion of external electrodes. This is shown in Fig. 4(b), 
where the y-component of the average spin accumulated 
on the dot is shown for several values of the parame- 
ters pi and pr- This figure clearly shows that the spin 
precession decreases when spin polarization of the leads 
becomes smaller. This is reasonable since lower spin po- 
larization creates smaller exchange field. 

When the spin polarization rate of the leads becomes 
smaller than 1, electric current can fiow also in the an- 
tiparallel configuration, contrary to the case shown in 
Fig. 2 for pi = pr = \. Consequently, the TMR effect also 
becomes smaller and remains finite for Lp = -k. In Fig. 5 
we show the TMR effect for different polarizations of the 
external electrodes, and for a particular noncollinear con- 
figuration [ip = 27r/3). The curves are symmetrical with 
respect to the bias reversal. The central peak corresponds 
to an enhanced magnetoresistance in the current block- 
ade regime. Such an enhancement of TMR in the block- 
ade regime, where sequential tunneling is suppressed, was 
also observed earlier for coUinear configurations [3. Al- 
though the sequential current is exponentially suppressed 
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FIG. 4: Bias dependence of the average spin component (Sy) 
for ip = 2n/3 and for indicated values of the Coulomb cor- 
relation parameter U (a) and the lead polarization (b). The 
other parameters are as in Fig.2. 
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FIG. 5: Bias dependence of the tunnel magnetoresistance in 
symmetrical junctions for indicated values of the lead polar- 
ization. The curves are plotted for the non-collinear configu- 
ration, ifi = 27r/3. The other parameters are as in Fig.2. 



in the blockade regime, the electrons still can flow due to 
higher order processes (like cotunneling). The other two 
broad maxima placed symmetrically on both sides of the 
central peak, occur between the two (lower and higher) 
threshold voltages. The following features of these two 
maxima are interesting to note. First, for small values 
of the polarization parameters, the central maximum is 
larger than the others. The situation changes with in- 
creasing polarization factors, and now the central peak 
becomes smaller for high spin polarizations of the leads. 
Second, the two maxima become strongly asymmetric for 
large values of the spin polarization factors, as clearly 
visible in Fig. 5. To account for this behavior one should 
take into account bias dependence of {Sy) from Fig. 4(b), 
which shows clearly the corelations between the spin pre- 
cession and height and shape of the TMR peaks with 
increasing polarization. 

In the situation studied above the dot level was empty 
at equilibrium. Qualitatively similar behavior of electric 

current, magnetoresistance, and average spin on the dot, 
has been found for the situation when the dot level is be- 



low the Fermi level, and the dot is occupied by a single 

electron at equilibrium. Exemplary current-voltage char- 
acteristics are shown in Fig. 6 (a), and the corresponding 
y-components of the average spin on the dot are shown 
in Fig. 6(b). The current- voltage curves display similar 
features as the curves shown in Fig. 2 (a), with characteris- 
tic negative differential conductance between the thresh- 
old voltages for noncoUinear configurations. The main 
difference in the bias dependence of {Sy) is a reversed 
sign of the peaks at the lower threshold voltage in com- 
parison to that in Fig. 3(c). 



B. Asymmetrical junctions 

Now the dot is separated from both electrodes by non- 
equivalent barriers, a ^ 1, and the electrodes are made 
of different ferromagnetic materials, pi i^Vr- For numer- 
ical calculations we assume pj=0.4, ^,.=1, and a = 0.1. 
More specifically, it is assumed that the right electrode 
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FIG. 6: Bias dependence of eleetric current (a), and the aver- 
age value of the y component of the dot spin (b), calculated 
for indicated values of the angle <p and ea = —0.1 eV. The 
other parameters are as in Fig.2. 



FIG. 7: Bias dependence of electric current (a) and magne- 
toresistance (b), calculated for indicated values of the angle 
ifi. The inset in (a) shows current between the threshold volt- 
ages for negative bias. The parameters assumed for numerical 
calculations are: Ed = 0.1 eV, U = 0.4 eV, Fo = 0.01 eV, 
pi = 0.4, pr = l,a = 0.1 and T = 100 K. 



is made of a half-metallic material with electrons be- 
ing totally spin-polarized at the Fermi level, whereas the 
factor Q = 0.1 indicates that on average electrons can 
tunnel much easier to (from) the left electrode than to 
(from) the right one. This asymmetry between the left 
and right electrodes and barriers gives rise to asymmetri- 
cal transport characteristics of the junction with respect 
to the bias reversal. The asymmetry is clearly visible 
in current-voltage characteristics and bias dependence of 
TMR, shown in Fig. 7(a) and Fig. 7(b), respectively. For 
positive bias (the right lead is the source electrode), the 
current and TMR curves are rather uniform above the 
first threshold voltage. The current flows there for arbi- 
trary value of the angle ip and thus TMR is significantly 
suppressed. The situation changes diametrically when 
the electric current flows in the opposite direction, i.e. 
when electrons tunnel through the dot from the left elec- 
trode to the right (half-metallic) one. Below the first 
threshold voltage sequential tunneling is exponentially 
suppressed and only the higher-order tunneling processes 



are possible. When the energy level ej, enters the tun- 
neling window, electric current starts to flow through 
the junction but this takes place only in a small voltage 
range in the vicinity of the first threshold voltage, where 
the characteristic resonant bump is observed. Above the 
bump, the current is suppressed by an electron residing 
on the dot. When Sd + U crosses the Fermi level of the 
source lead, the current increases again and finally sat- 
urates at a certain level. For positive bias, y > 0, the 
curves in Fig. 7(a) and Fig. 7(b) for diflferent values of the 
angle reflect a monotonous angular variation of the 
current and TMR, shown also explicitly in Fig. 8 for a 
particular value of the bias voltage. For negative bias 
one observes a more complex and interesting behavior of 
the transport characteristics. First, the above mentioned 
suppression of electric ciirrent between the two steps is 
now less pronounced, and the corresponding angular vari- 
ation of electric current and TMR is non- monotonous, as 
shown explicitly in Fig. 8 by the relevant curve. When 
the negative bias voltage surpasses the second threshold, 
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FIG. 8: Angular variation of electric current (a) and TMR (b) 
for indicated values of the bias voltage. The other parameters 
are as in Fig. 7. 



FIG. 9: Bias dependence of the average spin components; (Sz) 
(a), (Sx) (b) and (Sy) (c) for indicated values of the angle (p. 
The other parameters are as in Fig.6. 



the monotonous variation of electric current and TMR 
is restored. It is also interesting to note, that the TMR 
effect is enhanced in the voltage regions, where electric 
current is suppressed. This is the region below the first 
threshold voltage, and for negative bias also the region 
between the two threshold voltages. The latter one is par- 
ticularly interesting as the TMR may change there sign 
from positive to negative (see Fig. 7(b) and Fig. 8(b)). 

Suppression of electric current by an electron of a given 
spin orientation localized on the dot can be accounted for 
by analyzing spin accumulated on the dot when a steady 
state current flows through the system. This is illus- 
trated in Fig.9(a-c), where the part (c) shows the com- 
ponent induced by spin precession. The spin precession is 
particularly enhanced for negative bias between the two 
thresholds - exactly where electric current is suppressed. 
The enhancement is a consequence of relatively long time 
that electrons spend on the dot. 

The situation is different for the dot level lying under 
the Fermi level, i.e., when the dot is already occupied by 
one electron in equilibrium situation. The correspond- 



ing numerical results are shown in Fig. 10, where part (a) 
shows the current-voltage characteristics, and part (b) 
the average value of the perpendicular y-component of 
the average spin accumulated on the dot. The transport 
characteristics for collinear configurations were already 
accounted for in earlier publications [T^ , and their most 
interesting feature is the pronounced asymmetry with 
respect to the bias reversal (diode-like behavior). For 
noncoUinear configurations transport characteristics have 
new features, which are qualitatively similar to those 
found in the case of empty dot at equilibrium, like for 
instance nonmonotonous angular variation of electric cur- 
rent and TMR. 



VI. SUMMARY AND CONCLUSIONS 

Using the non-equilibrium Green function approach we 
have calculated electric current, average value of electron 
spin accumulated on the dot, and tunnel magnetoresis- 
tance due to rotation of the magnetic moments of exter- 
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nal electrodes. It has been shown that the average spin 
precesses by a certain angle around an effective exchange 
field arising from the interaction between the dot and 
electrodes. This precession leads to a nonzero value of 
{Sy), i.e., of the component normal to the plane deter- 
mined by magnetic moments of both electrodes. 



-'M ft6 



It has been also shown that the spin precession is en- 
hanced by Coulomb correlations and strong spin polar- 
ization of the leads. The spin precession is shown to 
exist also in the limit of vanishing Coulomb correlations 
on the dot. Moreover, the interplay of Coulomb correla- 
tions and effective exchange field may lead to a negative 
differential conductance in the voltage range between the 
two threshold voltages. It has been also shown, that the 
diode-like features of the system are partially suppressed 
when magnetic moments of the electrodes become non- 
collinear. 
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